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Abstract - The new classes of homogeneous cosmological models for the scalar fields are 
build in the context of Lyra's geometry. The different types of exact solution for the model 
are obtained by applying two procedures, viz the generating function method and the first 
. . - order formalism. 
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1 Introduction 



After the formulation of General Relativity (GR) by Einstein, many alternative geometric theories 

have been developed in order to explain gravitation phenomena. Inspired by geometrizing gravi- 

^•i ' tation, Weyl [T] proposed a more general theory in which both gravitation and electromagnetism 

are described geometrically. For a long time, Weyl's theory was not taken seriously due to non- 

O ' integrability of length of vector under parallel displacement. Later Lyra J2j suggested a modification 

of Riemannian geometry by introducing a gauge function which removes the non-integrability con- 

5-H ' dition of the length of a vector under parallel transport. This modified Riemannian geometry is 

^2j ■ known as Lyra's geometry. In contrast to Weyl's geometry, in Lyra's geometry the connection is 

metric preserving as in Riemannian geometry, and length transfers are integrable. It should be 

noted that Lyra introduced a gauge function into the structure-less manifold, as a result of which 

^ ' a displacement field arises naturally. Several authors (see, e.g. [3]-[Il]) have studied cosmology in 

Lyra's geometry. This alternating theory is of interest because it produces effects similar to those 

Q>«. ' produced in Einstein's theory 

(^ , Soleng [13] has pointed out that the cosmologies based on Lyra's manifold with constant gauge 

vector will either include a creation field and be equal to Hoyle's creation field cosmology [T4JI15J. 
f~^ , or contain a special vacuum field, which together with the gauge vector term, may be considered 

CNl ' as a cosmological term. Contrary to common assertion that the displacement vector field can play 

the role of the cosmological constant by itself, we want to emphasize that it can never play this 
role being alone, as seen from the equations presented below in this article. 

In general relativity Einstein succeeded in geometrising gravitation by identifying the metric 
^ tensor with the gravitational potentials. In the scalar tensor theory of Brans-Dicke on the other 

j^ , hand, the scalar field remains alien to the geometry. Lyra's geometry is more in keeping with the 

spirit of Einstein's principle of geometrisation since both the scalar and tensor fields have more or 
less intrinsic geometrical significance. Furthermore, the present theory predicts the same effects, 
within observational limits, as far as the classical Solar System tests are concerned, as well as tests 
based on the linearised form of the field equations [12] . 

To date, several authors have studied cosmology in Lyra's geometry with both a constant 
displacement field and a time-dependent one. For instance, in [12^ the displacement field is allowed 
to be time dependent, and the Friedmann-Robertson- Walker (FRW) models are derived in Lyra's 
manifold. Those models are free of the big-bang singularity and solve the entropy and horizon 
problems which beset the standard models based on Riemannian geometry. Recently, cosmological 
models in the frame work of Lyra's geometry in different contexts are investigated in several papers 
(see, e.g [Tl]-[21]). 

In last few decades there has been considerable interest in alternative theories of gravitation 
coursed by the investigations of inflation and, especially, late cosmological acceleration which is 
well proved in many papers |25j-|30j. In order to explain so unexpected behavior of our universe, 
one can modify the gravitational theory [3T]-[36|, or construct various field models of the so-called 
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dark energy (DE) which equation of state (EoS) satisfies w = p/p < — f/3. Presently, there is 
an uprise of interest in scalar fields in GR and alternative theories of gravitation in this context. 
Therefore, the study of cosmological scalar-field models in Lyra's geometry may be relevant for 
the cosmic acceleration models. 

Most studies in Lyra's cosmology involve a perfect fluid. Strangely, at least up to our knowledge, 
the case of scalar field in Lyra's cosmology was not studied properly. Here we would like to fill 
this gap. In this paper, we consider a scalar (quintessence or phantom) field and a tachyon field 
cosmological evolution in the context of Lyra's geometry. With motivation provided above, we 
have obtained exact solutions of Einstein's modified field equations for the spatially fiat Friedmann 
metric within the frame work of Lyra's geometry. For this purpose, we employ two methods, viz 
the generating function method and the first order formalism. 

2 Field equations 

The Einstein's field equations based on Lyra's manifold, as proposed in '3' and [T in normal gauge, 

may be written as 

1 3 3 

Rik - i^gikR + -;^4'i4>k - -7gik4>^(t>j = -n'^Tik, (i) 

where (/)i is the displacement vector, k^ = SttG and other symbols have their usual meanings in 
the Riemannian geometry. 

We assume a perfect fluid form for the energy-momentum tensor: 

Tik = {p + p)utUk -pgik, (2) 

and co-moving coordinates UiU* = 1, where Ui = (1,0,0,0). We also let (j)i be the time-like vector 

0, = (/3, 0,0,0), (3) 

where /3 = /3(t) is a function of time alone. The metric for FRW space-time is given by 

ds^ = dt^ - a^{t){dr^ + i'^{r)d^^), (4) 

where ^(r) = sinr, r, sinhr in accordance with a sign of the curvature k = +1,0,-1. For this 
metric together with 1^ and ([3]), the field equations ([T]) become 

3H' + ^- h^ = .^p, (5) 



2iJ + 3ij2 + — + -/32 = -K^p, (6) 



k 3 
a? 4 
where H = a/ a is the Hubble's parameter. 

Equations ^ and (H)) lead to the continuity equation as follows 



P+^/3/3 + 3^ 



-o + P+TT^/?" 



3 



0. (7) 



It is easy to find that the main equations of the model, i.e. ([S]), ^ and (O, can be presented in 
standard GR form, 

3fc 

3il' + — = K'pe//, (8) 



2H + 3H^ + ^ = ->i^Peff, (9) 



a^ 



and 



Peff + 3H 

by introducing two effective parameters: 



Peff + Peff 



(10) 



^3/32 3/32 



For the analysis of model, it is suitable to consider the EoS parameter w and the deceleration 
parameter q defined by 

Vpff 2 H a'd H , ^ 

In the absence of matter, that is when p = p = 0, the effective EoS is equal to w = +1, which 
corresponds to the so-called stiff fluid. That is why the displacement vector can never play the 
role of a cosmological term for which the EoS w = — 1 is required. 

To proceed further, we have to specify the type of scalar field. For the sake of simplicity, from 
now on we consider a fiat FRW cosmology: k — Q. 

3 Quintessence (phantom) Lyra's cosmology 

In this section, we consider a quintessence (or phantom) field as a source of gravity in Lyra's 
cosmology. Therefore, we have for the effective parameters (fTTj) as follows: 

Peff = ^P" + V{V) + g, Peff = |V^^ - V{^) + g, (13) 



where e = +1 represents quintessence while e = —1 refers to phantom field. In view of (|13|) . the 
set of basic equations ([8|),([9]) becomes 

3H^ = eip^ + 2V{^) + ^^^, (14) 

H = ^e^'-\p\ (15) 

where and in what follows we assume 47rG = 1. Even for the given potential V{tp) and displacement 
function /3(t), it is difficult to find exact solution for this model. However, a class of exact solutions 
can be obtained, say, in terms of the so-called generating function [371 131] ■ 



3.1 The generating function method 

Since /3(t) is an arbitrary function so far, we can transfer this arbitrariness into a new one, say 
coupling function f{(p{t)) as follows 

P' - p'i^)^'. (16) 

As a result, we have from p^ and (|16l) that 

Peff^l[fiv)+^]^^ + Vi^), (17) 

Peff = l[f{v^) + e]^^-Vi^). (18) 

Re-defying the scalar field as \P{^) + e|<^^ = i/"^, i.e. ip — ^ VlPiv) + ^1 ^^^ and putting 

e = sign[p{ip) -I- e], we get the following expressions for the effective parameters P^ : p^ff — 
s\'4>^ + U{4')i Peff = £\v ~ U{ip), where U{ip) — V{Lp{ip)). As seen in the effective parameters, 
the sign e of kinetic term can change (for e = — 1) from positive to negative and vice versa. 

Substituting ([TT]) and ([TSl) into Eqs. ([8|), ([9]), we have the following set of the main equations 
for our model: 

H' - l{\[fi^) + e]^' + Vi^)), (19) 

H = -[fi^) + e]^'. (20) 

In order to obtain exact solutions and following to the method of generating function [371 135] , we 
assume that F{ip) = (p. Applying the latter to the set of equations (rT9|) . pO)) . we can obtain the 
following general solution: 




Figure 1: The EoS parameter w, the deceleration parameter q and the Hubble parameter H in the 
case F{ip) = A versus x = At. Here, e = —1, wq = 5, f^ — 0.56, A = 4. 






and 



a{ip) = ooexp 






dip 



Moreover, as it follows from ([T2|) and (|22|) . the EoS parameter can be given by 
Let us consider two particular cases. 



(22) 
(23) 

(24) 
(25) 



3.1.1 The case F{ip) = A 

In order to obtain an explicit solution using this method, one has to specify the coupling function 
f{ip), or put some other additional condition. Here we assume that the coupling function can be 
represented by 

/(^) = 72/0 sin ^^. (26) 

With /q > 1/2 and e = —1, the kinetic term in p7)) . (IT51) can periodically change its signature, 
evolving from the phantom regime to the quintessence one and vice versa. From Eqs. ((2T|l - (|24p . 
we can find that 

ip — Af , 



H((^) = -A (/2 + £)(/.- ^ sinc^o'/' 

Wo 



P. 



a(</5) = ao exp I ^— — (^^ + —2^(1 - cos WQtp) |, 



and 



y{^)^% f (/o + e)^ - ^ sin 2c.o^^ ' ^' 
2 L Wo 



(/o +e)-/o coswovs 



(27) 
(28) 

(29) 
(30) 



The corresponding expression for the displacement function can be found from Eqs. p6p .(|26 p and 
((27|) as 

2^2, , . /woAi^ 



/3(t) = ^A/osm(^), 



(31) 




Figure 2: The EoS parameter w, the deceleration parameter q and the scale factor a in the case 
F{ip) — Xip versus x = Xt. Here, e = —1, (po = 1, /g = 1.12. 



According to ([72]) and (|28p . we have for the EoS parameter 

e + /o(l — coswoAt) 



w = — 1 



f 
(/o +e)Ai- —sin Wo At 

Wo 



2' 



(32) 



and the similar expression for the deceleration parameter q with 1 instead of the fraction 2/3. This 
solution is plotted in Fig. 1. It can be seen that the model expands in acceleration. At some 
points, where coswoAt — 1 — /q, the expansion is de Sitter {q = —1). The model periodically 
crosses the phantom divide in the both direction. 



3.1.2 The case F{ip) = Xip 

In order to obtain an explicit solution with some interesting feature in this case, we chose the 
coupling function f{(p) in the following form: 



For this choice, we have 



After integration, Eqs. (l2T])-([24l) gives 






ip = (pQ exp(Ai) 



H = -X 



"^0 



/o^At + /o^ln(po + e^exp(2Ai) 



a{t) = a, exp { - ^^ - f,Xt In^o - e^ (e^^* - l) }, 



(33) 

(34) 

(35) 
(36) 



where the constant of integration is chosen so that ao — a(0). Besides, we can find that the 
potential is given by 

V{^) = ^ f3(/2 In + i^2)2 _ ^^2 _ ^2l 



2 r^-'^"'^ 2 
After that, the Eos parameter can be found as follows 

/2 + e^2 cxp(2Ai) 



.. = -1 + - 



/2At + /2ln^o + e^exp(2At) 



(37) 



(38) 




Figure 3: The EoS parameter w, the deceleration parameter q, the scale factor a and the Hubble 
parameter H versus x — Xgt in the case e = +1, A = •\/3/2. Here, Xq = 3/2. 

The similar expression can be also obtained for the deceleration parameter. Some features of this 
solution is shown in Fig. 2. It can be seen that this model expands in acceleration as well. At the 
only point tg = A""'^ ln(/o/(/5o)j the expansion is de Sitter, that is q = —1, and the model crosses 
the phantom divide evolving from the quinessence sector to the phantom one. 



3.2 The superpotential method 

Other class of exact solutions can be obtained in terms of the so-called superpotential. This 
procedure was first performed with a single scalar field in 39 , and it was later re-opened as "the 
first order formalism" and extended on the case of two or more fields in [40]. The method of 
superpotential can be effectively applied to the quintom models as well (see, e.g., [41]). Keeping 
in mind the superpotential method, we now represent the geometrical field of displacement vector 
as a function of a new field a{t): 

P\t) = ^d?{t). (39) 



Now we get the new set of equations instead of p^ . (|15p : 



H 






Let us introduce the superpotential function W{ip, a) by the equation 

H^W{ip,a), 



(40) 
(41) 



(42) 



in which the Hubble parameter H{t), as a function of time, is presumably expressed in terms of 
fields (p{t), a{t). Substituting (|42|) into (|4T|) . one can obtain two first-order equations as follows: 



ip^ ~eW^, a = -Wa- 



(43) 



where W^ = dWjdkp, Wa = dW/da. The potential can be obtained from gO]), (|42]) and (gS]) as 



follows 



V{p) = -W^ 



:W,i 



1, 



:W' 



(44) 



2" 2' '^ 2' 

As the potential depends on ip alone, we have dV/da = 0. Taking into account (|43]l and (|44l) . it 
is easy to show that this is equivalent to 



dWg 
dt 



3WWa = 0. 



(45) 



The latter can be satisfied in many ways. Let us now consider one simple example of exact solution. 



We suppose that the superpotential is presented by 

W{^,a)^X{^)Y{a), 
that allows to rewrite p4|) in the following form: 

V{^) = \x^{^)[iY\a) - eY^{a)[^)" - Y'\a)_ 

We consider one example of the particular solution for this model based on 

X{ip) = Xoe-^'^, Y{a) = cosh{y/JIa), 
where ^ = 3 — eX^. So we have 

if, = -eW^ = eXXoe^^"^ cosh{y/JIa), 

a — ~Wa — ~^/JIXQe~ '^ sinh^y/JIa) , 
H{t) = Xoe^-^'^f*) cosh( V7Ia(t)). 
At the same time, we have the following expression for the potential (|m) : 



viv) 



^X^e-^V. 



(46) 
(47) 

(48) 

(49) 

(50) 
(51) 

(52) 



Combining and integrating Eqs. (|49p .(l5n i) . we can get 



exp(— Ai/?) — sinh(Y^a) 



(53) 



where the constant of integration is chosen to be equal to zero. Substituting the latter into (|50p . 
we have the following equation for a: 



= -y/JIXa sinh(y/Ia) 



M. 



(54) 



This equation can be integrated explicitly for several values of fi (or A and e) . 

As an example, we consider the case e = +1, A = \/3/2 or fi — 3/2. Integrating equation ([M)) 
and taking into account (j39p . we can obtain that 



a — \ — tanh 
3 



3Xot + 2 



f3' = 



32X§ 



[(3Xot + 2)2 - 4]2 



(55) 



where the constant of integration is chosen from the condition tanh{^3/2a(0)} = 1. With the 
help of ([5T|) . (|53p and ([55]) . it is possible to find that 



H{t) = 2Xo 



3Xot + 2 



{3Xot + 2) 



-, a{t)^ao {3Xot + 2) 



1/3 



and 



^ = 7!'" 



(3Xot + 2)2-4 



From Eqs. ([T^ and ([55)) . it can be found that 

4 



w 



1 



6 



(3Xot + 2)2' ^ 2 (3Xof + 2)2' 
Time evolution of this model is plotted in Fig. 3. 



(56) 



(57) 



(58) 



4 Lyra's cosmology of tachyon field 

In this section, we consider a tachyon field x ^-s a source of gravity in Lyra's cosmology. Substituting 
the well known expressions of tachyonic p and p into (|lip . we have 

Peff = 4^= + M!, p^^^ . -Vix)V^^+'-^- (59) 

As before, we put ff(t) = W[x{t),a{t)] and suppose the substitution (15^ . Due to ([5^ . the set of 
basic equations dS]), © becomes as follows 

ij = -il^-d^ (61) 

By inserting iJ = W{x,ct) into (1611) . one can obtain two first-order equations, 

where W;^ = dW/dx, Wa = dW/da. The potential is followed from (pO)) - (p^ in the form: 



^(^) = ^n ''' "'''' 



3VF2-M^2j _ 4^1/2 (g3) 



As this potential is independent on a, we have dV/da = 0. In view of (15^ . (p5|) . it is easy to prove 
that the latter is equivalent to equation (|^ . 

So we have a wide range of possibilities to solve the model equations assuming some certain 
dependence W{x,ct)- Instead, we can provide several classes of solution for the model evolving 
from some conditions on superpotential. Below, we show how it can be realized with the help of 
some ansatz for the superpotential. 

4.1 An example of exact solution for tachyon model 

One of the simplest ansatz for the superpotential may be written as follows: 

W^^X,W"'{x,a), W^^X2W"'{x,a), (64) 

where m, Ai, A2 are constants, and the equality W^^a — W^a,x i^ true. Therefore in view of (I45p 
and (|M|. we have 

H{t) = W{x{t),a{t)) = ^^^^, (65) 

where to is an integration constant. Hence, equations (j62l) become as follows 

2Ai(m/3)''"-^(t + to)'" . ^ (m/3)" 

^ 3(t + to)2™-2 _ A2(m/3)2"-2 ' " 2 (^ ^ ^^)„ . IDDj 



From the last equation in (|62|) and Eqs. (|64|) . ([65|) . it immediately follows that the displacement 
vector is given by 

/3(i)=A2^(y) (t + io)-™. (67) 

The EoS parameter and the deceleration parameter p^ are constant and defined by 

2 3 

«;==-l + — , <7 = -l + — . (68) 

777 777 

In order to distinguish between various DE models, Sahni et al. [42] proposed a cosmological 
diagnostic pair {r, s} called statefinder. The stateflnder test is a geometrical one based on the 
expansion of the scale factor a(t) near the present time to: 

a{t) = 1 + Ho{t - to) - \qoHlit - tof + \roHlit - tof + ..., 
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Figure 4: The points on {r, s} plan for different illustrative values of parameter m. The {1,0} 
point is related to the location of standard ACDM model. 

where a(io) = 1 and Ho,qo,ro are the present values of the Hubble parameters, deceleration 
parameter and the statefinder index r — a /aH^ respectively. The statefinder parameter s is the 
combination of r and q: s — (r — l)/3{q — 1/2). The important feature of statefinder is that the 
spatially flat ACDM has a fixed point {r,s} — {1,0}. Departure of a DE model from this fixed 
point is a good way of establishing the 'distance' of this model from flat ACDM. In terms of the 
Hubble parameter and its derivatives with respect to cosmic time the statefinder parameters of a 
flat FRW model are given by 

With the help of Eas. ([65)) . (|69|) . one can find that 

^ = l-- + 4' ^--' (70) 

for the model considered. The plot of {r, s} points corresponding to the different values of m is 
shown in Fig. 4. 

5 Conclusion 

In this paper, we have studied FRW cosmological models in normal gauge for Lyra's manifold with 
the quintessence (phantom) and tachyon scalar fields as the origin of gravity. We have built the 
new classes of FRW cosmological models of these scalar fields in the context of Lyra's geometry. 
The different types of exact solution for the model are obtained by applying two procedures: the 
generating function method and the first order formalism. We hope that the derived model is the 
next step in the development of Lyra's cosmology, and can be utilized to describe the dynamics of 
the evolution of the actual Universe. 
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